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Introduction 


In  this  report  we  investigate  the  extent  to  which  the 
vibrational  modes  of  a  system  of  combustion  effluents  consisting 
of  H2#  CO,  CC^,  and  HjO,  maintain  equilibrium  vibrational 
distributions  in  an  expanding  flow  field.  In  describing  the 
vibrational  relaxation,  only  the  first  vibrational  levels 
of  H2  and  CO  and  the  V2  (bending)  and  (asymmetric 
stretching)  modes  of  1^0  and  C02  are  included.  This  pre¬ 
cludes  any  quantitative  characterization  of  non-Boltzmann 

vibrational  distributions  such  as  are  expected  in  other 
(12) 

expansion  flows  '  but  a  discussion  of  their  possible 
effects  is  presented.  In  addition,  the  symmetric  stretch¬ 
ing  modes  of  the  polyatomic  species  are  not  explicitely 
included  primarily  due  to  a  lack  of  sufficient  rate  coeffi¬ 
cient  data  to  make  their  inclusion  meaningful. 

The  expansion  into  vacuum  of  flows  from  two  different 
nozzles  is  considered  in  this  report.  The  first  nozzle  is 
conical  and  no  strong  external  shocks  are  present.  The 
second  is  contoured,  and  a  strong  nozzle- induced  shock 
system  is  present  in  this  external  flow.  The  second  nozzle 
also  operates  at  considerably  higher  mass  flow  than  the 
first.  The  effect  of  vibrational  relaxation  on  the  flow 
enthalpy  has  not  been  included  in  these  calculations.  Since 
significant  vibrational  freezing  has  been  predicted  by  these 
calculations  the  results  must  be  viewed  as  qualitative 
particularly  in  the  outermost  streamlines.  However,  since 
many  of  the  important  vibrational  rate  coefficients  are 
not  well-known, particularly  the  vibration-vibration  exchange 
rate  coefficients  involving  polyatomic  species,  the  addi¬ 
tional  uncertainty  in  the  calculation  introduced  by  the 
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uncoupled  relaxation  treatment  is  probably  small. 

The  expanding  flow  fields  are  discussed  in  Section  I. 

In  Section  II,  the  details  of  the  relaxation  chemistry  included 
in  this  calculation  are  presented  together  with  the  rate 
coefficients.  In  Section  III,  the  relaxation  profiles  are 
presented  and  the  results  discussed.  The  major  conclusions 
of  this  analysis  are  presented  in  Section  IV. 
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Fluid  Mechanical  Calculations 


Because  the  temperatures  and  densities  within  the 
nozzles  examined  in  these  calculations  are  high,  we  expect 
significant  vibrational  lag  only  in  the  external  flow  fields. 
We  therefore  assume  that  the  flows  are  in  local  thermodynamic 
equilibrium  at  the  nozzle  exit,  and  use  the  nozzle  exit  as 
the  initial  station  in  our  calculations.  Upon  leaving  the 
nozzle,  the  gases  continue  to  expand  in  a  pattern  which 
depends  upon  the  flow  properties  and  the  exit  conditions, 
which  in  turn  depend  upon  the  nozzle  shape.  This  expansion 
continues  until  the  force  exerted  upon  th  ^anding  gases 
by  the  surrounding  medium  balance  those  ch  the  gas  exerts 
on  this  ambient  fluid.  Near  this  region  a  shock  structure 
develops?  in  order  to  avoid  the  complications  introduced 
by  this  interaction,  we  assume  that  the  expansion  is  into 
a  vacuum. 

We  also  assume  that  finite-rate  relaxation  processes 
have  a  negligible  effect  upon  the  thermodynamic  properties 
of  the  expanding  gases.  This  assumption  is  not  supported 
by  our  results.  However,  it  has  the  great  virtue  of 
allowing  us  to  decouple  the  relaxation  from  the  equations 
of  motion,  and  therefore  makes  the  total  calculation  much 
simpler.  In  addition,  it  does  not  affect  the  qualitative 
features  of  the  relaxation  calculation.  We  have  calculated 
the  flow  field  assuming  local  thermodynamic  equilibrium 
everywhere.  Vibrational  lag  would  result  in  a  somewhat 
lower  temperature  at  a  given  density,  and  in  a  slightly 
different  pattern  of  expansion,  and  therefore  an  enhanced 
tendency  to  "freeze"  the  vibrational  populations  at  a 
somewhat  earlier  time. 
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Two  separate  nozzle  flow  fields  have  been  examined  in 
these  calculations  (See  Figures  la  and  lb) .  Both  burn 
hydrocarbon  fuels  with  oxygen,  under  fuel-rich  conditions, 
and  the  major  product  species  are  l^#  1^0,  C0'  an<*  C02* 

The  first  nozzle  is  conical  (15°  half-angle) ;  the  flow 
rate  is  approximately  4  kgm/sec.,  with  an  assumed  exhaust 
composition  in  mole  per  cent  of  28.3%  ,  21.7%  1^0,  42.7% 

CO,  and  7.2%  CO2 •  We  assume  that  conditions  are  homogeneous 
and  source-like  at  the  nozzle  exit  plane,  where  the  tempera¬ 
ture  and  pressure  are  1550°K  and  0.68  atm.  The  expansion 

of  this  flow  is  calculated  with  the  aid  of  a  finite  differ- 

(3) 

ence  computer  program  allowing  for  variable  specific 
heats.  The  flow  field  is  shown  in  Figure  la  which  presents 
isotherm  in  the  external  field  and  also  shows  the  stream¬ 
lines  along  which  the  relaxation  calculations  were  performed. 
Temperature  histories  along  these  streamlines  are  shown  in 
Figure  2.  We  should  note  that  the  presence  of  a  wall 
boundary  layer  in  the  nozzle  would  have  a  significant 
effect  upon  the  flow  at  large  angles  as  well  as  on  the 
gas  properties  in  this  region.  We  have  not  included 
these  effects  in  our  calculations. 

The  second  nozzle  is  larger  than  the  first  and  incorpor¬ 
ates  a  number  of  additional  features.  The  nozzle  is  contoured, 
with  the  result  that  at  the  exit  plane  the  pressure  near 
the  wall  is  considerably  higher  than  at  the  nozzle  centerline. 
This  pressure  difference  causes  a  shock  system  to  form  in 
the  external  field.  (This  shock  system  is  produced  by  the 
nozzle,  and  is  separate  from  the  shock  system  which  develops 
as  the  expanding  flow  encounters  the  ambient  medium.  The 
nozzle  shock  system  is  present  even  in  the  flow  into  a 
vacuum.)  The  shock  created  as  the  wall  gases  expand  into 
the  central  gases  is  irregularly  reflected  from  the  axis, 
with  the  formation  of  a  Mach  disc.  We  also  allow  for  a 
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variation  in  flow  composition  normal  to  the  streamlines. 
(It  is  a  common  practice  in  nozzle  design  to  provide  a 
fuel-rich  layer  near  the  wall  to  reduce  thermal  loads  on 
the  structure.  Furthermore,  we  allow  for  the  presence 
of  an  annular  shroud  of  gases  of  different  compostion 
from  the  main  flow.  This  ducted  gas  confines  the  main 
flow,  and  prevents  its  expanding  to  larger  angles.  The 
total  main  flow  is  approximately  115  kgm/sec .  The  flow 
field  is  again  calculated  using  the  finite-difference 
procedure.  (The  calculation  is  discussed  in  detail  in 
Ref.  4) .  The  external  field  is  shown  in  Figure  16. 
Temperature  histories  along  these  streamlines  are  shown 
in  Figure  3.  Note  the  shock  structure  and  its  effect  on 
the  temperature  field  downstream  of  the  shocks;  we  shall 
see  that  this  has  an  effect  upon  relaxation  along  stream¬ 
lines  passing  through  the  shock  system. 


Vibrational  Relaxation  Chemistry 


In  this  section,  the  specific  excitation  processes  and 
rate  coefficients  used  in  these  relaxation  studies  are  presented 
The  extent  of  equilibration  of  the  vibrational  modes  with  the 
translational  temperature  as  given  in  Figures  2  and  3,  is 
calculated  for  the  species  CO,  H2,  C02  and  H20.  Only  the  first 
vibrational  level  is  included  for  the  diatomics  and  the 
symmetric  stretch  modes  of  the  triatomics  are  not  explicitely 
included  due  to  lack  of  sufficient  rate  coefficient  data  to 
make  their  inclusion  meaningful.  By  not  including  the  higher 
vibrational  levels  in  this  treatment,  the  possibility  of 
characterizing  the  non-Boltzmann  character  in  the  vibrational 
distributions  is  eliminated ^ .  However,  the  major  purpose  of 
these  calculations  is  to  estimate  the  extent  of  deviation  of 
the  vibrational  levels  from  equilibration  with  the  translational 
temperature  during  an  expansion  to  vacuum.  We  express  the 
results  in  the  form  of  effective  vibrational  temperatures 
through  the  equation 

(!)  n(l)  =  n (o)  exu{”E10/kTv1 

where  n(l)  and  n(o)  are  the  upper  state  and  ground  state  species 
densities,  respectively,  E1Q  is  the  vibrational  energy  level 
spacing  and  Ty  is  the  effective  vibrational  temperature.  The 
actual  vibrational  structure  and  the  appropriate  energy  spacing 
is  shown  in  Table  I. 

Vibrational  relaxation  of  exhaust  species  will  occur 
through  both  vibration- translation  (VT)  and  vibration-vibration 
(W)  energy  transfer  processes,  in  general,  w  processes  are 
rapid  compared  to  VT  processes  for  temperatures  below  about 
1000°K.  Thus,  some  vibrational  modes  with  relatively  slow  VT 
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TABLE  I 


COMBUSTION  SPECIES  AND  SPECTROSCOPIC  PARAMETERS 


SPECIES 


ENERGY  ABOVE  GROUND 


CO(m  = 

1) 

CO(v  = 

0) 

H2(v  = 

1) 

H2(v  = 

0) 

C02(v 

=  v3] 

C02(v 

=  V2) 

co2(v 

=  0) 

h20(v 

=  V3) 

h20(v 

=  V2) 

h20(v 

=  0) 

2143  cm"1 

4161  cm"1 

2349  cm"1 
667  cm"1 

3755  cm"1 
1595  cm"1 


TABLE  II 


RELAXATION  REACTIONS  AND  RATE  COEFFICIENTS 


REACTION 

VIBRATION- TRANSLATION  REACTIONS  RATE  COEFFICIENT 


1. 

CO(v=l)  +  h2  1  CO(v=0)  +  h2 

(1) 

^n(pr)  =  (68/T3)  -19.59 

2. 

CO ( v=l)  +  C02  «*  C0(v=0)  +  CO 2 

(2) 

*n(p r)  =  (241.16/T^)  -33.37 

3. 

CO(v=l)  +  H20  ?!  C0(v=0)  +H20 

(1) 

6.7~15exp[2.5~3T] 

4. 

H2(v=1)  +  H2  ?!  H2(v=0)  +  H2 

(3) 

Pt  =  3.9~10exp[100/T^l 

5. 

H2(v=1)  +  C02  ?!  H2(v=0)  +  CO 2 

(2) 

£n(pT)  =  (174.61/T3)  -21.49 

6. 

H2(v=1)  +  H20  ?!  H2(v=0)  +  H20 

(2) 

jen(pr)  =  (169.37/T^)  -21.37 

7. 

C02(v=v3)  +  H2  ?!  CO2(v=0)  +  H2 

(2) 

£n(pr)  =  (81.47/T3)  -19.85 

8. 

co2(v=v3)  +  C02  ♦*  C02(v=0)  +  C02 

(2) 

£n(pr)  =  (276. 27/T®)  -27.40 

9. 

C02(v=v3)  +  H20  ?!  C02(v=0)  +  H20 

(2) 

£n(pr)  =  (210.52/T%)  -24.40 

10. 

CO 2  ( v=v2)  +  Hj  jf  C02(v=0)  +  H2 

(1) 

3. 56-12exp[-8.04-4T] 

11. 

C02(v=v2)  +  C02  ?!  C02(v=0)  +  CO 2 

(1) 

5. 9-22T2 ' 77 

12. 

co2(v=v2)  +  h2o  ?!  co2(v=o)  +  h2o 

(1) 

3.95_11exp[-3.23"3T] 

13. 

h2°(v=V3)  +  H2  ?!  H20(v=0)  +  H2 

(2) 

jen(pT)  =  (147.7/T3)  -20.99 

14. 

H2°(v=v3)  +  C02  ?!  H2O(v=0)  +  CO 2 

(2) 

■tn (pr)  =  (393.48/T3)  -29.59 

15. 

H2°(v=v3)  +  H20  ?!  H20(v=0)  +  H20 

(2) 

£n(pT)  =  (330.28/T3)  -27.01 

16. 

H2°(v=v2)  +  H2  ?!  H20(v=0)  +  H2 

(2) 

£n(pr)  =  (47.16/T^)  -19.24 

17. 

h2o(v=v2)  +  co2  ?!  h2o(v=o)  +  co2 

(2) 

ln(pr)  =  (125.64/T^)  -21.99 

18. 

H20(v=v2)  +  H20  ?!  H20(v=0)  +  H20 

(1) 

8.74-12  +  9.36"9/T 

19. 

C02(v=v3)  +  H2  ?!  C02(v=V2)  +  H2 

(1) 

4-14 

20. 

C02(v=v3)  +  C02  ?!  C02(v=v2)  +  C02 

(1) 

3.58"14  +  9.43_17T 

21. 

C02(v=V3)  +  H20  ?!  C02(v=v2)  +  H20 

(1) 

5.47"13  -  3.03"16T 

22. 

H20(v=v3)  +  H2  ?!  H20(v=v2)  +  H2 

(2) 

jen(pT)  =  (70.66/T^)  -19.34 

23. 

H20(v=v3)  +  CO 2  ?!  H20(v=v2)  +  C02 

(2) 

Jtn(pr)  =  (188.24/T^)  -23.76 

24. 

H2°^v=v3^  +  H2°  ^  H2°^v=v2^  +  H2® 

(2) 

|n(pr)  =  (158/T^)  -22.53 
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TABLE  II  continued 


VIBRATION-  VIBRATION  REACT! ONS 


RATE  COEFFICIENT 


25.  H2(v=1)  +  CO(  v=0)  «♦  H2(v=0)  +  CO  ( v=l)  (4)  3 . 77~13T5/6exp[-88. 7/T^] 

26.  C02(v=v3)  +  CO(v=0)  «*CO2(v=0)  +  CO(v=1)  (4)  2~14T5/6exp[-36 . 55/T%] 

27.  H20(v=v3)  +  C0(v=0)  «♦  H2O(v=0)  +  CO(\^l)  (4)  3. 22-13T5/6exp[-136.1/T^] 

28.  H2(v=1)  +  H20(v=0)  M2(v=0)  +  H20(v=v3)  (4)  4.  ll“14T5//6exp[-31 . 5/T*] 

29.  C02(v=v3)  +  H20(v=0)  ^C02(v=v2)  +h2°(v=V2)  (4)  2.17“14T5/6exp[-23. 87/T^] 

30.  H20(v=v3)  +  C0(v=0)  £  H20(v=v2)  +  CO(v=l)  (4)  1. 34“15T5/6exp[-6.68/T3] 

31.  C02(v=v3)  +  H20(v=v2)  2  C02  ( v=0)  +  H20(v=v3)  (4)  1. 9“14T5/6exp[-30. 96/T%1 


13 


Pref erentially  relax  by  transfering  energy  to  molecules 
with  faster  VT  rates  through  W  coupling.  Although  only  a 
partial  set  of  vibrational  levels  have  been  included  in 
these  calculations,  an  attempt  has  been  made  to  include  all 
the  important  transfer  processes  which  affect  these  levels. 

^  relaxation  processes  included  is  given  in 

Table  II  together  with  the  adopted  rate  coefficient.  Since 
several  of  the  necessary  rate  coefficients  are  unknown,  they 
have  been  estimated  by  using  Millikan  and  White* s(6) 
correlation  for  VT  processes  and  from  the  modified  Rapp- 
Englander  Golden^  ^  model  for  W  processes. 

For  VT  rates,  the  anticipated  relaxation  time  values 
(pr)  have  been  converted  to  rate  coefficients  by  the 
equation ^ 


(2) 


1-e 


1 


where  T  is  the  translational  temperature.  No  major  uncertainty 
is  expected  from  the  use  of  this  conversion  equation  although 
the  estimated  values  of  p r  particularly  for  the  H20  molecule 
are  not  expected  to  be  very  reliable.  Fortunately,  some  new 
rate  coefficient  values  are  available  from  recent  molecular 
laser  systems  .  These  results  have  been  included  where 
appropriate. 

The  W  rates  have  been  either  taken  from  recent  measure¬ 
ments  of  interest  to  laser  system*9' 10)or  have  been  calculated 
from  the  following  rate  coefficient  expressions ^ 
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(3)  k  r,jj 


\,tol  =  1.508"1347/bl/  °/aa'i  Ea4^3  U?»  u»a*TB/i 


M-  “A 

H 


.  exp  {-1.32  (L,Ed)V3 
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/ 3 
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(4) 


vr » 

^3  0 


hdl 


=  8.88~la  p,  Ljf 


<T? 


mi* 


u; 


us  a 


{1  -  5 . 703-a  iiailEilli 

v  m  ' 


T 


where 

M-n  =  reduced  mass  in  molecular  units  ' 

La  =  exponential  parameter  in  A. 

O 

o a  =  hard  sphere  diameter  in  A. 

T  =  temperature  in  “K 
Ed  =  energy  defect  in  cirf1 

and 

^rn  —  matrix  element  for  r-»n  vibrational  transition. 

The  subscripts  on  these  W  rate  coefficients  refer  to  larger 

energy  defect  limit  (LEDL)  and  small  energy  defect  limit 

(SEDL)  and  suggest  the  range  of  application  for  each  equation. 

In  general,  when  the  difference  in  the  vibrational  quanta  being 

exchanged  is  larger  than  about  20  cm  \  equation  (3)  should 

be  used.  These  estimated  W  rate  coefficients  may  be  subject 

to  considerable  error  (although  probably  not  more  than  an 

order  of  magnitude  in  the  worst  case)  particularly  in  the 

temperature  dependence.  As  is  well  known,  the  temperature 

dependence  of  near  resonant  W  processes  may  become  negative 

at  temperatures  below  about  1000° K  due  to  long  range 

interactions.  This  has  been  confirmed  in  the  N„  -  CO 

2  2 
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(11) 

system.  Unfortunately  insufficient  data  are  available  to 
assess  this  possibility  for  other  systems  particularly 
involving  H20.  It  is  likely  therefore  that  the  temperature 
dependence  as  given  by  equations  (3)  and  (4)  will  not  be 
correct  for  very  near  resonant  W  processes  considered  here. 
However,  it  is  unlikely  that  the  qualitative  features  of 
the  relaxation  results  would  be  affected. 


16 


Relaxation  Results 


The  coupled  set  of  relaxation  equations  with  the 
appropriate  rate  coefficients  as  given  in  Table  II  have  been 
integrated  along  the  various  streamlines  shown^ ^Figures 
2  and  3,  using  a  Runge-Kutta-Merson  procedure  '  .  The 

translational  temperature  and  density  were  treated  as 
known  functions  of  time  for  the  integration. 

The  Runge-Kutta-Merson  integration  procedure  used 


features  a  variable  time  step  consistent  with  equation 
stability  and  also  contains  an  algebraic  algorithm  for 
integrating  those  relaxation  equations  which  are  sufficiently^ 
close  to  steady  state  to  make  the  set  of  equations  "stiff 


Representative  results  for  the  relaxation  calculation 
are  given  in  Figures  4  through  8  for  the  conical  nozzle 
flow  and  in  Figures  9  through  13  for  the  contoured  nozzle 
flow.  AS  expected,  the  low  lying  bending  modes  of  both  C02 
and  H-Q  show  significant  vibrational  freezing  only  in  the 
outermost  streamlines  for  which  the  expansion  is  most  rapid. 
The  higher  lying  vibrational  modes,  particularly  of  the 
diatomics,  show  freezing  behavior  for  almost  all  streamlines 
considered.  For  those  cases  where  the  relaxation  behavior 
is  calculated  through  the  Mach  disc,  i.  e.  Figure  9  , 
the  flow  remains  in  equilibrium  until  some  distance  downstream 
from  the  Mach  disc  due  primarily  to  the  high  densities  in 


this  region  of  the  flow. 
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Figure  4:  Vibrational  temperature  history  along  the  center  stream¬ 
line  of  the  conical  nozzle  flow  field.  Hie  assumed 
composition  is  28.3%  H2,  21.7%  H20,  42.7%  CO  and  7.2%  CO, 


If 


TIME  (SEC) 


Figure  5:  Vibrational  temperature  history  along  a  streamline  in  the 
conical  nozzle  flow  field  for  which  m  =  0.0394  mt.  The 
assumed  composition  is  28.3%  H_,  21.7%  H_0,  42.7%  CO  and 
7 . 2%  CO_ .  2 
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Figure  9:  Vibrational  temperature  history  along  a  streamline  just 
outside  the  Mach  disc  (m  =  0.878  mt>  for  the  contoured 
nozzle  flow  field.  The  assumed  composition  is  10.3%  Hj 
40.4<&  H-O.  25.2%  CO  and  24.1%  CO- . 
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Figure  11:  Vibrational  temperature  history  along  a  streamline  in 

the  contoured  nozzle  flow  field  for  which  m  =  0.663  mfc. 
The  assumed  composition  is  10.3%  H_,  40.4%  H_0,  25.2%  CO 
and  24.1%  C0„.  22 
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Figure  13:  Vibrational  temperature  history  along  a  wall  streamline 
of  the  contoured  nozzle  flow  field  for  which  m  =  0.066  m 
The  assumed  composition  is  23.5%  26.1%  1^0,  39.4%  CO 

and  11.0%  C0„. 
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Discussion 


The  major  conclusions  to  be  drawn  from  these  relaxation 

calculations  is  that  the  bending  modes  of  K20  and  CO2  are 

in  equilibrium  with  the  translational  temperature  in 

essentially  all  regions  of  the  expanding  flow  except  the 

outermost  streamlines  where  the  expansion  is  most  rapid. 

This  result  is  not  expected  to  be  affected  by  improved 

rate  coefficients  as  enough  information  is  available  on 

the  kinetics  of  the  bending  modes  to  make  their  modeling 

reliable.  The  stretching  modes  of  HjO  and  C02  and  the 

diatomic  modes  are  largely  out  of  equilibrium  and,  in  fact, 

the  diatomic  modes  are  normally  frozen  out  at  or  near  the 

nozzle  exit  conditions.  Since  the  coupling  rates  between 

(9) 

the  v3  and  v2  modes  of  C02  are  moderately  well  known' 

(i.e.  reactions  (19)  -  (21)),  the  calculations  involving 
the  v3  mode  are  expected  to  be  fairly  reliable.  The  comparable 
rate  coefficients  involving  the  H20  molecule  (i.e.  reactions 

(22)  -  (24))  are  not  known.  However,  the  predicted  rate 
coefficient  for  reaction  (22)  is  of  the  same  order  as  the 
analogous  reaction  involving  C02 (reaction  (19)).  This 
appears  reasonable  since  vibration-rotation  energy  transfer 
tends  to  make  the  H20  molecule  relax  faster  than  CC>2  but  the 
larger  quantum  being  transferred  in  the  H20  case  tends  to 
counter  balance  this  effect.  The  estimates  of  reactions 

(23)  and  (24)  are  not  expected  to  be  reliable  to  even  an 
order  of  magnitude.  This  was  not  considered  a  serious 
problem  since  the  concentration  of  H2»  H20,  and  CO^  are 

all  comparable,  and  therefore  reaction  (22)  dominates  as  over 

either  (23)  or  (24) .  No  experimental  rates  are  available 

on  the  W  processes  but  from  analogy  to  known  W  rate 

(9)  (13) 

coefficients  in  C02  -  N2  and  N2  -  CO  systems  ,  equa¬ 
tions  (3)  and  (4)  are  expected  to  give  rather  reliable 
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values  at  any  given  temperature.  As  indicated  earlier, 
however,  the  temperature  dependence  for  the  near  resonant 
W  processes  (i.e.  reactions  (29)  and  (30)  and  perhaps 
even  for  reactions  (26)  and  (31) )  may  be  negative  at  low 
temperatures.  This  would  favor  improved  coupling  between  the 
Vj  and  modes  and  would  suggest  somewhat  lower  vibrational 
temperatures  for  the  mode  than  shown  in  Figures  (4) 
through  (13) .  Due  to  the  rapid  decay  in  the  translational 
temperature  and  density  for  these  vacuum  expansions,  it  is 
unlikely  that  improved  calculations  would  alter  the  qualita¬ 
tive  features  of  the  figures. 

Finally,  we  note  that  vibrational  modes  remain  in 
equilibrium  with  translation  in  passing  through  the  normal 
shock  (Mach  disc)  formed  by  the  contoured  nozzle.  However, 
in  the  re-expansion  downstream  from  the  Mach  disc,  signifi¬ 
cant  departure  from  equilibrium  is  found.  This  departure 
will  be  greater  for  smaller  nozzles,  where  the  expansion  is 
more  rapid. 

In  conclusion,  our  calculations  show  appreciable 

vibrational  freezing  of  combustion  effluents  particularly 

in  the  outermost  streamlines  of  both  the  conical  and  contoured 

nozzle  flow  fields.  This  conclusion  is  not  expected  to 

bo  altered  by  improved  rate  coefficients  although  quantiative 

results  will  vary  as  both  W  and  VT  processes  are  more 

accurately  determined  for  those  reactions  involving  HjO. 

Due  to  the  predicted  high  vibrational  populations  and  the 

low  translational  temperatures,  non-Boltzmann  vibrational 

distributions  should  be  evident  as  a  result  of  anharmonic 

(1.2) 

W  pumping 
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